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ABSTRACT 

Spicules are ubiquitous, fast moving jets observed off-limb in chromospheric spectral lines. Combining the 
recently-launched Interface Region Imaging Spectrograph with the Solar Dynamics Observatory and Hinode , 
we have a unique opportunity to study spicules simultaneously in multiple passbands and from a seeing free 
environment. This makes it possible to study their thermal evolution over a large range of temperatures. A recent 
study showed that spicules appear in several chromospheric and transition region spectral lines, suggesting that 
spicules continue their evolution in hotter passbands after they fade from Ca II H. In this follow-up paper we 
answer some of the questions that were raised in the introductory study. In addition, we study spicules off-limb 
in C II 1330 A for the first time. We find that Ca II H spicules are more similar to Mg II 2976 A spicules than 
initially reported. For a sample of 54 spicules, we find that 44% of Si IV 1400 A spicules are brighter toward the 
top; 56 % of the spicules show an increase in Si IV emission when the Ca II H component fades. We find several 
examples of spicules that fade from passbands other than Ca II H, and we observe that if a spicule fades from a 
passband, it also generally fades from the passbands with lower formation temperatures. We discuss what these 
new, multi-spectral results mean for the classification of type I and II spicules. 


1. INTRODUCTION 

Spicules have a long and rich scientific history, but their 
relatively small widths, barely resolved even in modern tele¬ 
scopes, combined with a complex dynamical behavior make 
spicules a source of lively scientific debate on their origin and 
role in the solar atmosphere. Traditionally, spicules are ob¬ 
served above the limb in chromospheric spectral lines such as 
Ca II H 3968 A and H<x In the traditional picture, spicules are 
seen to move outward before falling back toward the limb a 
few minutes later. Typical lifetimes observed are 2-6 min with 
typical velocities of 15-40 km s -1 . Athay & Holzer (1982) 
suggest that the release of gravitational energy from falling 
spicules could be the source of heating in the atmosphere and 
Pneuman & Kopp (1978) indicated that significant spicule 
material could reach coronal temperatures. Withbroe (1983) 
studied EUV emission from the downfalling spicule material 
and indicated that release of gravitational energy is not the 
primary source of heating in the upper chromosphere. See 
Beckers (1968, 1972) for reviews of older spicule observations 
and Sterling (2000) for a review on the theoretical side of early 
spicule models. 

With the emergence of modern space borne telescopes and 
telescopes utilizing advanced image post-processing, a wealth 
of new discoveries has been made. De Pontieu et al. (2007) 
used Ca II H data from the Hinode satellite to discover a 
different class of spicules, most common in quiet Sun regions 
and coronal hole regions. These so-called type II spicules are 
predominantly seen to move outward from the limb before 
fading from the observed passband. Their lifetimes are smaller 
than what was traditionally observed for spicules, usually less 
than 2 minutes. The few constraints on the formation of type 
II spicules, together with their possible vital contribution to 
the energy budget (due to their sheer number density) have 
placed type II spicules at the center of many debates, both 
regarding their very existence and their formation (De Pontieu 
et al. 2007; Judge et al. 2011; Zhang et al. 2012; Pereira et al. 
2012, 2013; Lipartito et al. 2014; Pereira et al. 2014). 


One possible explanation for the fading of Ca II spicules 
is heating. De Pontieu et al. (2011) show that the disk coun¬ 
terpart of type II spicules can be connected to brightenings in 
the hotter passbands of the Atmospheric Imaging Assembly 
(AIA, Lemen et al. 2012) onboard the Solar Dynamics Obser¬ 
vatory, which might make the new class of spicules vital in 
the energy balance of the outer solar atmosphere. Pereira et al. 
(2014, hereafter Paper I) use co-observations from three differ¬ 
ent spacecrafts and find that type II spicules have companion 
components in the Mg II 2976 A and Si IV 1400 A filtergrams 
of the Interface Region Imaging Spectrograph (IRIS, De Pon¬ 
tieu et al. 2014b) in addition to He II 304 A from AIA. The 
spicules continue to evolve after fading from the Ca II H pass- 
band, strongly suggesting that spicules are the site of heating 
to at least transition region temperatures. 

Numerical simulations have yet to shed any significant light 
on this predominantly observational feature, with only one 
notable exception (Martinez-Sykora et al. 2011). However, 
numerical experiments will be invaluable in determining in 
detail the mechanism for formation and heating of spicules. 
Providing observational constraints will help verification and 
interpretation of simulation efforts and is one of the goals of 
this paper. We tackle this issue by looking for answers to the 
following questions: How common are spicules that disappear 
from Ca II H and brighten up in Si IV at the same time? Are 
Mg II spicules really that different from Ca II spicules? What 
fraction of Si IV spicules are brighter near the top? How do the 
decelerations of type II spicules measured in hotter passbands 
compare with those reported for type I spicules? 

2. ANALYSIS 
2.1. The observations 

In the remainder of this paper the Ca II H 3968 A, 
Mg II 2976 A, Si IV 1400 A, C II 1330 A and He II 304 A 
filtergrams are referred to as Ca II, Mg II, Si IV, C II and He II, 
respectively. 

We use coordinated observations from instruments onboard 
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Fig. 1.— Limb cutouts of the two datasets analyzed. The red dashed line is at the same position in all images of a given dataset. The photospheric limb is nearly 
parallel to the x-axis at y ~ 2 Mm. The thick vertical black line seen in the IRIS slit-jaw images is the IRIS slit. In the Ca II images we can see the border between 
the two CCD sensors as a vertical line at x ~ 20 Mm because we are pushing the data and looking at slight differences in the dark current correction, which is not 
perfect. All images except He II are radially filtered (see text). 


three different space borne telescopes: the Broadband Filter 
Imager (BFI) on the Solar Optical Telescope (SOT) onboard 
Hinode (Tsuneta et al. 2008; Suematsu et al. 2008; Kosugi 
et al. 2007), AIA and IRIS. From Hinodel BFI we analyze 
filtergrams centered on the Ca II H spectral line. The filter 
bandwidth is 0.3 nm, the exposure time is 1.8 s and the pixel 
size is O'/ll pix -1 . From IRIS we analyze the following slit- 
jaw (SJI) filtergrams: SJI 2796, centered on the line core of 
Mg II, SJI 1330, dominated by the C II lines, and SJI 1400, 
dominated by the Si IV lines. The far-UV passbands, C II and 
Si IV, have a bandwidth of 40 A, and the near-UV Mg II has 
4 A band pass. All IRIS slit-jaw images have 8 s exposure time 
and a pixel size of O'/17 pix -1 . From AIA we use the 304 A 
channel which is dominated by He II. The exposure time is 
2.9 s and the pixel size is 0'/6 pix -1 and the full width at half 
maximum (FWHM) of the filter is 12.7 A. 

Assuming equilibrium conditions we expect the Ca II and 


TABLE 1 
IRIS DATA SETS. 


Starting time 

(, x , y ) Coord. 

At 

Duration 

Filtergrams 

2014-02-21T1L24 

(7", 

-987") 

19 s 

90 min 

Mg II, Si IV 

2014-02-19T16:28 

(6", 

-965") 

23 s 

67 min 

C II, Si IV 


Note. — At is the cadence of the Mg II or C II images (for Si IV the 
cadence was always 19 s). The duration covers the overlapping period between 
Hinode and IRIS. The first dataset is the same as in Paper I. 


Mg II passbands to be dominated by radiation from plasma 
at 10,000 to 15,000 K, C II from 15,000 to 25,000 K, Si IV 
from 60,000 to 80,000 K and He II from 75,000 to 100,000 K 
(Golding et al. 2014; Carlsson & Leenaarts 2012; De Pontieu 
et al. 2014b). 

We use two different datasets with a different set of IRIS 
filtergrams available for each, see details in Table 1 . The first 
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dataset is the same as in Paper I. At the time of the observations 
no coronal hole was present on the South Pole so the region 
is characterized by quiet Sun conditions. The cadences of 
BFI/SOT and AIA were 5 s and 12 s, respectively. 

We made extensive use of CRISPEX (Vissers & Rouppe van 
der Voort 2012) when studying the observations. 

2.2. Aligning the datasets 

The Hinode and AIA data are interpolated to the IRIS plate 
scale. The IRIS data were rotated by 0.6 degrees relative to 
the center of the slit-jaw images to compensate for the slightly 
different orientation between IRIS and AIA. 

The images observed with the various passbands are differ¬ 
ent to a varying degree. The alignment of the data from the dif¬ 
ferent spacecrafts poses a formidable challenge because there 
is no common filtergram or one that is sufficiently similar to be 
used for unambiguous automatic alignment. Because we are 
analyzing limb observations, we expect that filtergrams formed 
at different heights to have a spatial offset in the limb-ward 
direction. We chose the 1600 A channel from AIA to serve as 
a common base for the alignment. The height of formation of 
the channels is different in magnetic and non-magnetic regions. 
On-disk in non-magnetic regions, the 1600 A, Si IV and C II 
channels are dominated by continuum radiation originating 
from the upper photosphere, and the Ca II and Mg II passbands 
are formed at comparable heights (M. Carlsson, private com¬ 
munication, see also Martmez-Sykora et al. 2015). In magnetic 
regions, however, in particular the Si IV and C II channels, 
and to some extent the 1600 A channel, are dominated by line 
emission from transition region lines, and formation heights 
are much less confined. The procedure is then to align IRIS 
and Hinode data to each other indirectly through the 1600 A 
channel in non-magnetic regions, because there we expect the 
channels to be formed at roughly the same height. A precise 
alignment between the IRIS channels is made possible by the 
fiducial mark on the slit (De Pontieu et al. 2014b) and the He II 
passband is aligned to the 1600 channel to within a precision 
of 1 or 2 AIA pixels by internal alignment with the use of level 
1.5 AIA data (Lemen et al. 2012). 

After various experiments we settled on the following ap¬ 
proach: The Si IV to 1600 alignment was done separately 
for the x and y direction. For the y direction we used cross¬ 
correlation on a subsection of the disk field of view (FOV), 
approximately 470 Mm 2 , strictly avoiding network elements. 
For the x direction we used cross correlation on the entire 
disk-part of the IRIS data excluding the slit in the slit-jaw 
images. To reduce projection effects we put the sub-FOV as 
far from the limb as the data allowed. Cross-correlation on 
two other, non-magnetic, sub-FOVs was done to estimate the 
uncertainty of the alignment. We performed the x direction 
cross-correlation on the entire disk-part, including network, 
elements, as it gave overall better visual alignment compared 
to only using the sub-FOV. When using the entire disk-part 
for the y direction, the computed shift in y direction would be 
dominated by the magnetic network regions, for which TR line 
emission results in an ill-determined offset in the y direction. 
This effect was not visible in the x direction, likely because 
the network-based shifts canceled out because the network 
features are predominantly oriented along the y axis, not the x 
axis. 

The Ca II to 1600 alignment was performed manually. Shifts 
were applied every 285 s, and the shift of the intervening 
frames were computed from linear interpolation. We estimate 


the uncertainty to be within 077, approximately the same as the 
AIA resolution. We are not able to achieve sub-pixel accuracy 
because the images in the channels are not similar enough to 
compare large-scale patterns. We estimate that the uncertainty 
is the same for the Si IV to 1600 alignment. We remark that 
after numerous experiments we had to accept a relatively large 
uncertainty in the final alignment of the datasets. For on- 
disk multi telescope coordinated observations, we have been 
able to achieve much more accurate alignment to sub-pixel 
precision (e.g., De Pontieu et al. (2014a); Martmez-Sykora 
et al. (2015); Rouppe van der Voort et al. (2015)). For future 
coordinated limb observations we suggest the inclusion of a 
pure photospheric channel such as the IRIS 2832 SJI and the 
G-band filtergram on BFI/SOT. 

We used a temporal grid with 19 s cadence and for each time 
selected the images closest in time in the different passbands. 
19 s was chosen because that was the cadence of the Mg II and 
Si IV data. For the time of each image we used the middle 
time of the exposure. 

2.3. Description of the FOV in the different filtergrams 

Figure 1 shows an overview of the limb part of the dataset 
with a dashed red line as a visual aid. For both datasets the 
red line is on top of a spicule. All other thin (less than 5 Mm 
wide) features seen sticking out are considered to be spicules. 
In the first dataset the red line is on top of a wide structure, 
wider than most spicules, and in the Ca II and Si IV images 
we can see hints of substructure, multiple “threads”, inside 
the large structure. All the “threads” move synchronously in 
time. Skogsrud et al. (2014) studied similar “threads” in higher 
resolution Ha data from the Swedish 1-m Solar Telescope and 
this may be the manifestation of those threads in the Ca II and 
Si IV passbands. 

Our general impression from working with these filtergrams 
is that the Mg II images share characteristics with both the 
Si IV and Ca II images, making them an excellent diagnostic 
for detecting spicules. The large scale structures are similar to 
the structures in the Si IV image, while many of the small scale 
features are similar to Ca II structures, and the combination of 
Mg II and Ca II made it easier to find more spicules in the first 
dataset compared to the second dataset, which did not include 
Mg II. Many spicules are barely observed in Si IV, for example 
in the region between x = [20,40] Mm in the top example of 
Figure 1 where many clear spicules are visible in Ca II and 
Mg II. 

The He II image also appears to show many spicules, but the 
spatial resolution of AIA and potentially the optical depth are 
a clearly a limiting factor in resolving small scale structures. 
We also note that the average height at which the lowest point 
of spicules can be discerned is greatest in the optically thick 
lines of He II and Mg II, while considerably lower in Ca II and 
Si IV. 

The spicule traced by the red dashed line in the first dataset 
is clearly visible in all images except Si IV, where it appears 
very faint. In contrast, the spicule sticking out to the right from 
the dashed line appears strongest in Si IV. 

In the second dataset, which contains C II but not Mg II, we 
observe that many spicules in C II, appear faint and close to 
the noise level in the data. There is little spicule substructure 
visible in the C II slit-jaw images. 

2.4. Radial filters 

Ca II H spicules, as observed by Hinode, are extremely faint 
and barely stand out in the original filtergrams. Radial density 


4 


filters are then employed to enhance their visibility (e.g., De 
Pontieu et al. 2007; Zhang et al. 2012; Pereira et al. 2012; 
Paper I). These filters are built by dividing each filtergram by 
(/(r)), the normalized mean spatial and temporal intensity as 
a function of r, the distance from the limb. This greatly in¬ 
creases the contrast between the spicules and their background. 
However, (/(r)) is mostly noise for large values of r, and 
dividing the filtergrams by such values becomes meaningless. 
Therefore, when building the radial filter (I(r)) is scaled so 
that for large values of r it tends to a constant value a that is 
higher than the background noise. The choice of a is largely 
subjective and determined by visual inspection, so that spicules 
are enhanced but artefacts and noise are minimized. In Paper I 
and earlier work, a = 0.135 was used. 

Spicules in the IRIS filtergrams are not as faint as in Ca II H, 
but radial filtering is still important to improve their visibility. 
In addition, as noted in Paper I, spicule trajectories are easier 
to follow in IRIS than in Ca II H because they do not fade as 
much. The analysis of space-time diagrams in Paper I showed 
that after fading, Ca II spicules still showed an extremely faint 
outline that seemed to match the evolution of spicules in other 
filters. In this work we investigated this point further and 
found that this residual spicule signal in Ca II filtergrams can 
be enhanced by using a more aggressive radial filter. Using 
a = 0.018 we get noisier images farther from the limb, but 
spicules are seen to extend farther out, to about the same 
heights as the Mg II or Si IV spicules. The point of the 
radial filter is an important one because, as discussed in earlier 
literature, the height and lifetime measurements of spicules are 
critically dependent on their appearance as seen in filtergrams. 

The two radial filters are compared in Figure 2. The 
very faint spicules seen sticking out from about (x,y) = 
(14,15) Mm toward the upper left of the image are not identi¬ 
fiable in the non-filtered intensity, and are barely discernible 
with the radial filter of Paper I (see slight bump in intensity in 
the lower panel). However, they are seen more clearly with 
the new filter. By comparing with the intensity along the black 
dotted column (serving as an indicator of background inten¬ 
sity level, where no spicules exist at such heights), one can 
appreciate that the spicule signal is very low, comparable to 
the noise level (black dotted line). The temporal evolution of 
this spicule is shown in Figure 2 of Paper I. Making a time 
series figure equivalent to their Figure 2 using our improved 
radial filter, we can trace the Ca II spicule component for much 
longer: up to 456 s with the improved filter, compared to 248 s 
with the old filter. 

We applied the improved radial filter to the Ca II images. 
For the IRIS images there was no advantage in using this more 
aggressive filter, so we used the same radial filter as in Paper I. 
For the AIA He II images no radial filter was applied. 

2.5. Space-time plots 

The primary tool we use for analyzing the temporal evolu¬ 
tion of spicules is space-time plots, xt -plots, with x being the 
distance along a slit and t time. To obtain an xt -plot of the 
spicules, a slit was manually drawn on the spicules for each 
time step for a 10 minute duration. This allows for the tracking 
of the transverse motion of spicules. The slits were drawn on 
the Ca II and Mg II images. Due to the alignment uncertainty, 
we need to take into account that the position of the features 
can be slightly shifted between the passbands. To accommo¬ 
date these small shifts, a slit width of l'/5 was used. For each 
height in the extended slit, the maximum intensity value was 
computed from the radial filtered images. We chose the maxi- 
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Fig. 2.— Comparison of radial filters. Top panel: image with the radial filter 
used in Paper I (inverted color table). Middle panel: image with the radial 
filter used throughout this work. Bottom panel: comparison of radial filters 
applied above. The dashed lines show the radial filter intensities (/(r)) for 
the filters of Paper I (gray) and this work (black). The solid gray/black lines 
show the intensities along the vertical lines of the same color in the panels 
above (gray: the radial filter of Paper I, black: the radial filter of this work), 
and the solid blue line shows the original non-filtered intensities along the 
same vertical line. The dotted black line is the intensity along the dotted black 
line in the middle panel. The arrows indicate the top of a spicule. 


mum value instead of the average value to avoid the artifacts 
which then can occur. For example, if we move the artificial 
slit from one frame to the next to track a narrow spicule, and 
the spicule is perfectly centered in the first slit, but is slightly 
offset in the next, the average value will be shifted to a lower 
value because more of the dark background is included in the 
slit. It was verified that using a smaller or wider slit had little 
effect on the xt- plots. 

The statistics for the spicules were obtained from the xt- 
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plots in all spectral passbands. The lifetime was defined to be 
the time difference between the onset and end of the visible 
spicule, i.e. the lifetime only measures the visible part of 
the spicule in the xt- plots. The lifetime of a spicule differs 
in the different passbands and the lifetimes we measured are 
a lower estimate of the actual lifetime of the physical event: 
we can only trace isolated spicules above the “spicule forest” 
and most likely part of the lifetime of most spicules is spent 
hidden between other spicules. The minimum height at which 
the spicules are visible varies between the passbands, and 
between spicules due to interfering neighboring spicules. For 
example, in Ca II the lowest positions of the spicules can often 
be seen around 4 Mm above the photospheric limb, while for 
Mg II the lowest position is often seen at around 8 Mm above 
the limb. The way we detected spicules, by searching in all 
passbands, may give a bias toward the longer lived spicules. 
Our cadence of 19 s is also not sufficient to detect spicules with 
lifetimes of less than about a minute. In Pereira et al. (2012), 
the lifetimes were determined based on filtergram inspection. 
It is possible that by using xt -plots instead, with the knowledge 
of the spicule position after it fades in Ca II, that a bias toward 
longer lifetimes is introduced because the faint trails left after 
the spicule fades might not be associated with the spicule in 
xy -snapshots at the time. 

Our length measurement follows the spicule from the top, at 
time of maximum extent, down toward the limb to the point 
where the line-of-sight superposition, “spicule forest”, ob¬ 
scures the spicule. From that point the shortest line to the 
photospheric limb is added. If the spicule has an angle to the 
limb normal, or the actual footpoint is below the limb due to 
the projection in the observations, we are underestimating its 
length. 

We measured the peak ascent velocity by fitting a straight 
line to the steepest part of the xt- plot. We applied an edge 
enhancement algorithm beforehand to the Ca II xt -plot to 
make the manual extraction easier for the spicules that did not 
have a clearly defined upper point. Because there is a large 
uncertainty in the manual approach, we did the measurements 
twice to estimate the consistency of our results. The average 
difference between the two measurements was -2 km s -1 with 
a standard deviation of 21 km s _1 . Figure 3 illustrates how 
the velocities were extracted with two different good fits with 
about 17 km s _1 difference between them. 

3. RESULTS 

We measured a total of 54 spicules in the two datasets, 35 in 
the first and 19 in the second. The large difference in number of 
spicules detected is because the second dataset was of shorter 
duration, but also because Mg II was not observed, which 
made cross-channel identification harder (see Section 2.1). 

In Figure 4 we show the time evolution of two spicules, one 
from each dataset. The first spicule appears to have significant 
substructure and is seen rising and falling in all passbands 
except Ca II, where it fades significantly. The slit is centered 
on the left-most of the two substructure “threads”. The length 
of the spicule is comparable in all passbands. For the second 
spicule the length is also comparable in the passbands. The 
spicules fade in the downfall phase in the Ca II, C II and Si IV 
passbands. The C II component appears very faint. 

In Paper I, three spicules are shown in detail and all are 
clearly detected in Ca II, Mg II, Si IV and He II. From the 
spicules we measured in the same dataset, we find that 23 
out of 35 spicules (66%) have a clear component in all four 
passbands. Of the 12 that do not have a clear component in 



Fig. 3.— Illustration of how the peak ascent velocity was extracted from the 
Ca II xt -plots. The upper panel shows an edge-enhanced version of the Ca II 
image in the bottom panel. The two dashed lines in each panel both represent 
an equally good estimate of the velocity (values printed in the legend). 


all passbands, the Si IV component is unclear in 10 and the 
He II component in 2. By unclear we mean that the relatively 
thin structure of the spicule seen in the other passbands is not 
identifiable, instead we observe emission from a much larger 
structure, possibly obscuring the spicule. Of the 54 Ca II 
spicules we analyze we find that 41 (76%) fade rapidly, usually 
in less then 38 s (i.e., in two timesteps). From those that fade, 
22 (54%) show a faint trail in the xt -plots at the location of 
the spicule in passbands other than Ca II. These faint leftover 
traces are not easily connected to the earlier spicule phase if 
we were not aided by knowledge of the position of the spicule 
in the other passbands, see also Paper I. Of the Mg II spicules, 
11 out of 35 (31%) are seen to fade, 3 out of 36 (8%) Si IV 
spicules fade, 1 out of 19 (5%) C II spicules fade and only 1 
out of 54 (2%) He II spicules are seen to fade. Mg II spicules 
typically fade over a height range between 3-10 Mm in less 
than 38 s. We observe that if a spicule fades from a passband, it 
generally also fades from the passbands with a lower formation 
temperature. For example, a spicule can fade from the Mg II 
passband and still be visible in the Si IV passband, but not the 
other way around. 

In Figure 5 we show xt -plots for six different spicules se¬ 
lected to illustrate the typical behavior seen in our data. For 
each spicule we show a snapshot, in different passbands, from 
the time given by the vertical dotted red line in the xt -plots on 
the same row. The bottom three spicules are from the second 
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Fig. 4.— Time series of two spicules, one from each dataset. The spatial cutouts are fixed in space and the artificial slit used for the construction of the 
^-diagrams is between the two red dashed lines. The corresponding xt -plot to the top spicule is #1 in Figure 5 and #4 for the bottom spicule. Note that all images 
except He II are radially filtered. 


dataset, containing C II. For spicules #l-#5 there is a visible 
component in all passbands, and for #6 only the Ca II, C II and 
Si IV components are clearly visible. For all the Ca II spicules 
we observe that the final phase is significantly weaker com¬ 
pared to the onset phase, #3 and #4 are very clear examples. 
In spicule #3 we observe emission from only the top of the 
spicule in Ca II during downward phase after it has faded. The 
path traced by the top matches the parabola seen in the other 
passbands. For the Si IV components we see that #l-#3 and #5 
show an increase in emission in the final phase and in #2 and 
#5 the Si IV component is brighter toward the top. We find that 


24 of the 54 (44%) Si IV spicules show increased emission in 
the top part of the spicule. 

In Paper I it is reported that, for some examples, the Si IV 
emission is increased when the Ca II spicule component fades. 
We find that in 30 out of the 54 spicules (56%) the Si IV chan¬ 
nel displays increased emission when the Ca II components 
fade. Spicules #l-#3 and #5 in Figure 5 are clear examples 
of this behavior. In #1 and #2 we see that the Si IV compo¬ 
nent brightens up when the Ca II components fade, and that 
the Si IV xt -plot is a natural (parabolic) continuation of the 
Ca II xf-diagram. By looking at the evolution of intensities we 
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Fig. 5.— xt-diagrams of six spicules from the two datasets. The first column is a spatial cutout of the FOV at the time marked with a red dashed line in the 
corresponding xt -plots on the same row. The slit regionis between the two red dashed lines. The black and white stripes visible in Ca II for spicules #3 and #6 are 
due to missing data. Note that all images except He II are radially filtered. Time series of spicule #1 and #4 are shown in Figure 4. Animations of this figure is 
available in the online material. (Movies here: movie-1, movie-2, movie-3, movie-4, movie-5, movie-6 
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Peak velocity [km/s] 


Fig. 6.— Comparison of our results with previously reported statistics for 
quiet Sun linear spicules (QSL), type II spicules, in Ca II. In the bottom panel 
we included reported velocities from active region parabola spicules (ARP), 
type I spicules. Note that the ARP velocities were computed from the first 
derivative of the parabolic fit and the other velocities are from fitting a straight 
line to the xt -plots. 

estimate the delay between the channels to be about 3 minutes. 
In #3 and #5 we observe that more of the middle part of the 
Si IV spicule brightens up below the top point when the Ca II 
component disappears. Table 2 summarizes selected results. 

It is of interest to compare the statistics of our Ca II data with 
the existing numbers reported for type II spicules in Ca II. In 
Figure 6 we compare our measured spicule properties with the 



Deceleration [km/s^2] 

Fig. 7.— Comparison of our measured Mg II decelerations with previ¬ 
ously reported decelerations for active region parabola spicules (ARP), type I 
spicules, in Ca II. 


TABLE 2 

SPICULE EVOLUTION STATISTICS 


Element 

# spies 

% that fade 

% show parabolic trf-path 

Call 

54 

76 

33 

Mg II 

35 

31 

82 

C II 

19 

5 

64 

Si IV 

54 

8 

68 

He II 

54 

2 

62 


results of Pereira et al. (2012) for quiet Sun spicules in Ca II 
data. The figure shows histograms of the maximum height, 
peak ascent velocity and lifetime. The lengths we measure 
are between 8-20 Mm with 12.5 Mm average for Ca II and 
13 Mm average for Mg II. The difference in results between 
the Hinode , IRIS and AIA passbands are small. In the lifetime 
plot we can see that the Ca II lifetimes are slightly shifted 
toward shorter duration. This is what we expect based on the 
observation that more spicules fade from Ca II compared to 
the other passbands. The Ca II average lifetime is 5 min, while 
for the remaining passbands it is about 6.5 min. The longest 
lifetime found is 10.5 min. The average peak velocity we 
measure is 53 km s -1 and the maximum velocity we find is 
108 km s -1 . We find that the lengths and lifetimes we measure 
are significantly longer than previously reported. The veloc¬ 
ities match well for quiet Sun spicules, and are significantly 
higher than reported type I spicule velocities. 

In Figure 7 we compare our measured Mg II decelerations 
with previously reported numbers for active region parabola 
spicules, identified as type I spicules. The decelerations match 
well with reported type I decelerations. We determined the 
deceleration from fitting a parabola to the Mg II xt- plot. The 
reason for choosing Mg II is that the Ca II and Si IV xt -plots 
often show an incomplete parabola, while the Mg II shows 
a complete parabola. We performed the measurements on a 
subset of 22 Mg II spicules that showed clear and complete 
parabolas in the xt -plot for its entire duration. 

4. DISCUSSION 

We studied the detailed evolution of 54 spicules in the 
Ca II H, Mg II 2976 A, Si IV 1400 A, C II 1330 A and 
He II 304 A filtergrams. We find that a typical spicule, in 
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quiet Sun regions, has a visible component in all filtergrams. 
The Ca II component typically fades away and the evolution 
continues in the other passbands. Most Ca II spicules that fade 
leave a faint “trace” in space-time plots. We find that 44% of 
the Si IV spicules are brighter toward the top, while 56 % of 
the spicules show an increase in Si IV emission when the Ca II 
component fades. We interpret this behavior as the effect of 
heating in the spicule, which is consistent with earlier work 
(De Pontieu et al. 2007; De Pontieu et al. 2011; Pereira et al. 
2012; Paper I; Rouppe van der Voort et al. 2015). 

In Paper I it was found that Mg II and Si IV spicules continue 
to rise for 4-8 Mm after the Ca II spicule component fades 
and that the Mg II and Si IV spicule lifetimes were longer by 
several minutes. In this study, using a more aggressive radial 
filter on the Ca II data, we find that Ca II spicules are more 
similar to Mg II spicules than it appeared in Paper I. We are 
pushing the Ca II data to the noise limit and exploiting the 
access to Mg II data to detect very faint spicule signal that 
was otherwise indiscernible. We observe much more similar 
heights and lifetimes for Ca II and Mg II spicules. This is 
what we expect because the elements have roughly the same 
formation temperature. The small difference in lifetimes we 
find between Ca II and Mg II are perhaps due to the different 
opacities of the spectral lines, Mg is about 18 times more 
abundant than Ca (Asplund et al. 2009). However, because the 
intensity ratio of Ca II and Mg II is not constant as function of 
time there must be other effects that also play a role. 

We find several examples of spicules that fade from pass- 
bands other than Ca II, and we note that if a spicule fades from 
a passband, it also generally fades from the passbands with 
lower average formation temperature. We find more spicules 
that fade from the cooler passbands than from the hotter pass- 
bands. 

When comparing our results to Ca II quiet Sun spicules of 
Pereira et al. (2012) we find good agreement in peak ascent 
velocity. When comparing our results to type I spicules, the 
active region parabola spicules of Pereira et al. (2012), we find 
similar decelerations for our Mg II spicules as for Ca II type I 
spicules, but we find much higher peak velocities in the ascent 
phase than in Ca II type I spicules. 

It is important to discuss what these new multi wavelength 
results on spicules mean for our understanding of the difference 
between type I and II spicules. Now that we find clear parabolic 
paths in the Mg II and higher temperature passbands and at 
least traces of parabolic paths in part of the Ca II spicules, the 
question arises whether distinction between types of spicules is 
still justified. The spicules we study have higher peak velocity 
compared to type I spicule, and they clearly undergo heating, 
which is not compatible with type I spicules. The longer 
lifetimes and heights we can extract from the data do not alter 
the fundamental aspect of type II spicules, that they fade in 
the different passbands and particularly in the cooler (or lower 
opacity) passbands. We observe that the Ca II spicules fade 
significantly and that the downward phase is primarily only 
visible in the IRIS and AIA data. The parabolic paths in the 
Ca II data only become apparent after aggressive filtering and 
it pushes the data to the limits fundamentally set by noise. 

It is now clear that the spicules we study are subject of 
heating. If we assume that ionization equilibrium is valid in 
these dynamic events, then the standard estimated temperature 


ranges for the different passbands would be correct, and we 
would conclude that we always observe spicules to consist of 
a mix of chromospheric and transition region temperatures. 
However, the dynamic nature of these jets indicates that non¬ 
equilibrium ionization may well play an important role, which 
would change our interpretation of what occurs in spicules. 
Regardless of the ionization state, our observations support a 
scenario in which heating is occurring, often toward the top 
of the spicule. We measure typical timescales between peak 
intensity in Ca II and Si IV to be about 3 minutes. This can 
be interpreted as a typical time scale for the heating process, 
although it is possible that non-equilibrium ionization effects 
could change this value. 

In our space-time plots we generally observe parabolic mo¬ 
tions. The Ca II xt-diagrams often show only the early, more 
linear, phase of the parabola. Parabolic paths can be the result 
of driving by shock waves (Hansteen et al. 2006; Heggland 
et al. 2007, 2011). However, it is difficult to envision a sce¬ 
nario in which magneto-acoustic shocks (of the type seen in 
Hansteen et al. (2006)) could lead to observations of heating 
like in our observations. In shocks the heating is occurring 
at the shock front as it passes through the material, and the 
aftermath of the shock is responsible for the parabolic motion. 
This does not seem to be compatible with our observations. 
The Mg II parabolas are generally the most complete, unlike 
in Ca II which often only shows the beginning, and in Si IV 
which is sometimes weak in the beginning and shows increased 
emission in the later stages. One can speculate whether this is 
an optical depth effect: possibly the spicule is optically thick 
in Mg II, but optically thin in Ca II and Si IV. Heating of the 
spicule, not necessarily uniformly across and along the whole 
spicule structure (e.g., in a subset of threads), might result 
in fading from the optically thinner Ca II while there is still 
sufficient opacity in Mg II in cooler parts of the spicule. Only 
the hotter part of the spicular plasma can then be seen in Si IV. 
This observational puzzle puts strong constraints on numerical 
models. 

The increased emission we see in the Si IV xt- plots, often 
occurring in the downward phase of the spicules, is compatible 
with previous on-disk observations of an apparent correlation 
between red-shifted H-alpha spicular plasma and C IV emis¬ 
sion (de Wijn & De Pontieu 2006). It also makes one wonder 
what role spicules play in the pervasive redshifts observed in 
transition region lines. This warrants further investigation. 


IRIS is a NASA Small Explorer mission developed and op¬ 
erated by LMSAL with mission operations executed at NASA 
ARC and major contributions to downlink communications 
funded by the NSC (Norway). Hinode is a Japanese mission 
developed by ISAS/JAXA, with the NAOJ as domestic partner 
and NASA and STFC (UK) as international partners. It is 
operated in cooperation with ESA and NSC (Norway). The re¬ 
search leading to these results has received funding from both 
the Research Council of Norway and the European Research 
Council under the European Union’s Seventh Framework Pro¬ 
gramme (FP7/2007-2013) / ERC grant agreement nr. 291058. 
B.D.P. is supported by NASA contract NNG09FA40C (IRIS), 
and NASA grants NNX11AN98G and NNM12AB40P. 
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